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ABSTRACT 

In the majority of the literature on plasma shock waves, electrons play the role of “ghost particles,” since their contribution to mass 
and momentum flows is negligible, and they have been treated as only taking care of the electric plasma neutrality. In some more 
recent papers, however, electrons play a new important role in the shock dynamics and thermodynamics, especially at the solar-wind 
termination shock. They react on the shock electric field in a very specific way, leading to suprathermal nonequilibrium distributions 
of the downstream electrons, which can be represented by a kappa distribution function. In this paper, we discuss why this anticipated 
hot electron population has not been seen by the plasma detectors of the Voyager spacecraft downstream of the solar-wind termination 
shock. We show that hot nonequilibrium electrons induce a strong negative electric charge-up of any spacecraft cruising through 
this downstream plasma environment. This charge reduces electron fluxes at the spacecraft detectors to nondetectable intensities. 
Furthermore, we show that the Debye length A 1 ^ grows to values of about A^/Ad - 10 6 compared to the classical value Tp in this 
hot-electron environment. This unusual condition allows for the propagation of a certain type of electrostatic plasma waves that, at 
very large wavelengths, allow us to determine the effective temperature of the suprathermal electrons directly by means of the phase 
velocity of these waves. At moderate wavelengths, the electron-acoustic dispersion relation leads to nonpropagating oscillations with 
the ion-plasma frequency cu p , instead of the traditional electron plasma frequency. 
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1. Introduction 

The majority of the plasma-physics literature on shocks es¬ 
sentially considers general flux-conservation requirements only, 
lea ding to the well-known Rankine-Hugoniot relations (e.g. . 


seelSerrin ll959tlHudsonlll970t iBaumiohann & Treumannf l 996; 


Cfombosil 1998 : [Diver 2001 ). These relations, however, do not 

explicitly formulate the internal microphysical processes that 
generate internal entropy during the conversion from the up¬ 
stream regime into the downstream plasma. The missing physics 
in this description evidently leads to the well-known phe¬ 
nomenon that the set of Rankine-Hugoniot relations is a math¬ 
ematically unclosed system of equations. Therefore, these rela¬ 
tions can only provide unequivocal solutions if additional phys¬ 
ical relations are added to the system, such as the assumption 
of an adiabatic reaction of the plasma ions during their compres¬ 
sion into the higher-density regime on the downstream side (e.g., 
lErkaev et alJl2000h . 

The solar-wind termination shock is a particular example of 
a plasma shock for which microphysical effects play an impor¬ 
tant role. According to recent studies, pick-up ions have a cru¬ 
cial influence on overall shock physics at the solar-wind termina¬ 
tion shock. They are a thermodynamically important additional 
plasma component, since they extract a significant fraction of the 
upstream kinetic ene rgy in the form of t h ermal energy at the ter¬ 
mination shock (see lDecker et al.l [20081. Z ank et all (1201 Oh and 


iFahr & SiewerJ (2 0071 201(1 . 201 1) have studied kinetic features 
of this multicomponent shock transition and found relations be¬ 
tween upstream and downstream ion distribution functions that 
are different for solar-wind protons and pick-up protons. Al¬ 
though these studies discuss the required overadiabatic reaction 
of pick-up ions, a satisfying explanation of all plasma properties 
observed by Voyager-2 ( Richanisonetal .120081) is still lacking as 
demonstrated bv lchalov & Fahr (2013). The latter authors show 
that assuming a significant difference in the behavior of solar- 
wind electrons compared to protons, namely as an independent 
plasma fluid, leads to an explanation of most of the observed 
plasma data presented by iRichardson et al . d2008i) in a satisfy¬ 
ing manner. 

To achieve this result in their parameterized study, 
IChalov & Fahrl (120131) had to include preferential heating of 
the solar-wind electrons during the shock passage by a fac¬ 
tor of about ten stronger than the proton heating. This type of 
electron heating at the potential jump of fast-mod e shocks had 
been r ealize d earlier bvlLer ov & M angenevl (Il984l) . iTokar et ail 


(1986J), and ISchwartz et a" ( 19881) . and the phenomenon of 


shock-heated electrons also appears in plas ma-s hock simulations 
when electrons are treated kinetically (see Lembege et alJl2003l 
12004)) . In these cases, the plasma electrons demagnetize due to 
two-stream and viscous interactions and attain downstream-to- 
upstream temperature ratios of 50 and more. 
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iLerov et all (Il982l) and iGoodrich & Scudded (1 19841) follow 
a different approach. In their treatments, the plasma electrons 
carry drifts perpendicular to the shock normal (z-direction) dif¬ 
ferent from the ion motion as a reaction to the shock-electric 
field. These drifts establish an electric current j ± , which is re¬ 
sponsible for the change of the surface parallel magnetic field 
B\\ in the form 4 nj, /c = dBy/dz. To achieve the same consis¬ 
tency, i Fahr et aT| (j2012j ) describe the conditions of the upstream 
and downstream plasma in the bulk frame systems with a frozen- 
in magnetic field. In this framework, the Liouville-Vlasov theo¬ 
rem describes all relevant downstream plasma quantities as an 
instantaneous kinetic reaction in the velocity distribution func¬ 
tion during the transition from upstream to downstream. The 
excessive electron heating is then the result of the mass- and 
charge-specific reactions to the electric shock ramp, as shown 
in the semikinetic m odel s of the multifluid term ination shock 
bv IFahr et al.l (1201 2l) and IFahr & Siewertl (120131) . According to 
these studies, electrons enter the downstream side as a strongly 
heated plasma fluid with negligible mass density that dominates 
the downstream plasma pressure. 

In this paper, we demonstrate why the Voyager-1/-2 space¬ 
craft did not detect these th eoretically suggested hot electrons 
(see Ikichardso n et al.l 120 08) when they penetrated into the he¬ 
liosheath plasma. For the purpose of clarification, we analyze 
the downstream plasma conditions in more detail under which 
the detection of preferentially heated electrons would have to 
take place. 


2. Theoretical description of downstream electrons 

In the following section, we shall start from a theoretical descrip¬ 
tion of solar-wind electrons expected downstream of the termi¬ 
nation shock iFahr & Siewertll2013l) . We treat them as a separate 
plasma species, which reacts in a very specific manner to the 
electric-field structure connected with the shock before adapting 
to the downstream plasma bulk frame. In the shock-at-rest sys¬ 
tem, the shock electric potential ramp decelerates the upstream 
protons from the upstream bulk velocity U i to the downstream 
bulk velocity (A.p, which is comparable to the center-of-mass 
flow U* — (/o,p( 1 + 5 y/m/M), where 5 is the shock compression 
ratio and m, M denote the masses of electrons and protons, re¬ 
spectively. The downstream magnetic field is frozen-in into the 
center-of-mass flow, and all plasma components are eventually 
comoving with the center-of-mass flow iChashei & Fahrll2013l) . 
The electrons, on the other hand, react in a completely different 
way to this electric potential. First, they attain a strong “over¬ 
shoot” velocity f/ 2 , e which then relaxes rapidly to the center- 
of-mass bulk velocity (7* enforced by the frozen-in magnetic 
field. During this relaxation process, the plasma generates ran¬ 
domized thermal velocity components through the action of 
the two-stream instability or the Buneman instability as well 
as by pitch-ang le scattering (see IChashei & Fahrl 1201 3l i2014t 
IFahr et alfzoTi . 

Under the assumptions of an instantaneous reaction of the 
electrons to the electric potential and randomization of the over¬ 
shoot energy by the Buneman instability and pitch-angle scatter¬ 
ing to a n isotropic distribution in th e downstream bulk frame, we 
obtain (IFahr & SiewertfTo 13ll2015l) the following expression for 
the electron pressure P 2 , e on the downstream-side of the shock: 


B 2,e 


M s 1 - 1 U\ 
m 5 cr 

l,e 


|A(a) sin 2 a + B(a) cos 2 o-J Pp p . 
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( 1 ) 


The indices p and e denote proton- or electron-relevant quan¬ 
tities, and the indices 1 and 2 denote upstream and downstream 
quantities, respectively. The parameter U denotes the bulk veloc¬ 
ity, and 5 is the shock compression ratio. The parameter c \ >e is the 
average electron thermal velocity on the upstream side, where 
we assume that solar-wind electrons and protons have equal tem¬ 
peratures. We denote the magnetic tilt angle between the shock 
normal and the upstre am mag netic field as a . The functions A{a) 
and B(a) are given in lFahr & Siewerl)(l2013l) . In the limit of van¬ 
ishing thermal pressures and dominating magnetic pressures, the 
self-consistent compression ratio s turns out to be s — 1 (see 
IFahr & Siewerjf20l3l Eq. 18). Instead of treating the calculated 
velocity moment P 2 , e of the distribution function, we focus on 
the distribution function / 2 , e itself to derive the expected electron 
particle fluxes g 2 ,e as the relevant observable for the Voyage r-1/- 
2 instrumentation. As motivated in IFahr & Sicwcrt d'2013h . we 
assume that the downstream nonequilibrium distribution func¬ 
tion /i.e is a general kappa-distribution. This convenient choice 
represents the transition from a thermal core to a suprathermal 
tail distribution. The kappa-distribution is given by 


/e(v) = 


r(* e + 1 ) 


^3/2 3/2q 3 r(*e - 1/2) 


K e ®~ 


-te+1) 


( 2 ) 


where n e is the electron number density, 0 e is the velocity width 
of the thermal core, and K e is the specific electron kappa parame¬ 
ter with a range of 3/2 < A^e, < The symbol I = r(x) denotes 
the Gamma function of the argument x. 

In the next step, we determine the adequate value of K 2 >e 
for electrons downstream of the shock associated with a pres¬ 
sure given by Eq. (Q}. For this purpose, we determine the as¬ 
sociated pressure P 2 ,c,k (he-, the pressure resulting as the sec- 
ond velocity moment o f the above distribution, Eq. 0, see 
IHeerikhuisen et al1l2008i) that is equal to the pressure given by 
Eq. (D (i.e., the electron pressure found in the multifluid ap¬ 
proach bv IFahr & Siewertl ( 2013 ) ). We obtain the following re¬ 
lation for P 2 ,e,K- 


Pi,c,k(k) — 3/2.c _ Go ( 


*2,e 


K2,e - 3/2 


M s 2 - 1 u\ 

111 5 c? 

l,e 


[ A{a ) sin 2 a + B(a ) cos 2 ar]Pi p . 


(3) 


As shown in lFahr & Siewertl (12013i) . we can define the factor If 
by 


9 3KT i p 3Pi p 

0 2 e = n-— = n— 


m 


n\ p m 


(4) 


where Pi p is the upstream solar-wind proton pressure. The factor 
If describes the change of thermal core velocities from upstream 
to downstream. We quantify this factor later in Section^ Using 
an upstream proton temperature of Pi p = 2 x 10 4 K and 11 ~ 1, 
we obtain an average energy for the thermal core electrons of 
<e e ,c) = (l/2>n0 2 e = KT ,, p = 1.72 eV. 

From the above relation Eq. (0, we first obtain: 


"2,e n 3 f l.p K2,e 
2 «l,p *2,e-3/2 

Mr-1 t/? 9 , 

=- tt~\M a) sin” a + B(a) cos 2 a]P\ p , (5) 

m 5 cr 

l,e 
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which further simplifies to 


*2,e 2 1 M S 1 - 1 U\ 

K 2e - 3/2 3 Tim s 2 c 2 c " 


( 6 ) 


in the region behind the shock, where U 2 , e ~ Ui$, using the short 
notation A(a) = A(a ) sin 2 a + B(a ) cos 2 a. 

Relating the thermal velocity ci >e of the upstream electrons to 
the thermal velocity ci iP of the upstream protons through T i iP = 
T\ fi , we obtain 


K2,e 

K2,e - 3/2 


2 s 2 - 1 

3n s 2 


U\ 

~T~M a)- 

cf „ 


(7) 


We assume that the upstream solar-wind Mach number of the 
protons (i.e., /q -p = f/i/ci iP ) is of order 8, which then leads to 


3 

2i^8 2 A(a)-3/2' 


( 8 ) 


The termi nation-shock compressio n ratio observed by Voyager-2 
is s ^ 2.5 ((Richardson et alj[2008l) . which leads to 


3 54A(a)/n 

~ 2 54A(a)/n - 3/2' 


In the case of a perpendicular shock (i.e., a - 90°), A(7r/2) = 
A(7 t/2 ) = s leading to 


3 135/n 

2 135/11 - 3/2’ 


( 10 ) 


Assuming a value of n = 1, we obtain the result k 2c = 1.517. 
This kappa index k 2 # for the shocked downstream solar-wind 
electrons characterizes a highly suprathermal electron spectrum 
with a power-law nearly falling off as v~ 5 as shown in Eq. (0. 
Consequently, we can write for the resulting distribution func¬ 
tion of the downstream electrons. 


/ 2 ,e(v) 


«2,e 


r(K 2 ,e + 1) 


n 2 ' 2 ^ 2 ® 2 r(* 2 ,e 

2,e 2,e 


1 / 2 ) 


1 + 


K2.C& 2 „ 


~(*2,e + l) 


( 11 ) 


with /C 2 ,e — 1.517. 

The distribution function in Eq. CU is easily transformed 
into the spectral electron flux g 2 ,e(v) = 47rv 3 /2, e (v). We normalize 
velocities as x = v/© 2 ,e and find 


g 2,eW 


4H 2 ,e r(/C 2 , e +1) , 

n 1 ' 2 K 3 ' 2 IX*2,e - 1/2) X 

2,e 



— (^2,e + l) 


( 12 ) 


In Figure [I] we show these spectral electron fluxes for different 
indices k 2>e . 


3. Electric equilibrium potential 

In the following, we calculate the electric equilibrium poten¬ 
tial up to which a spacecraft charges up when entering the 
heliosheath under the assumption that both electrons and pro¬ 
tons have kappa-distribut ions and that emission processes are 
negligible in the plasma. (Fahr & Siewerd d2Q 1 3t> show that all 
ions (solar-wind as well as pick-up ions) treated as one fluid 
can be characterized as one joint kappa-function with a joint 
kappa-index k 2 i — 2, depending on the pick-up ion abundance 
downstream of the shock (see Figure 2 of IFahr & Siewertl201 3i) . 



Fig. 1. Normalized spectral electron fluxes downstream of the solar- 
wind termination shock according to Eq. GD. We show the kappa- 
distributed fluxes for different values of k 2 ,c- The special case /C 2 , e = 
1.517 is our result for the perpendicular solar-wind termination shock 
according to Eq. ( 1 101) . 


Therefore, we describe electrons with Eq. (ITU and protons with 
the following distribution downstream of the shock: 


/ 2 ,i(v) 


«2.i 


n> 2 ,i +1) 


^/ 2 ^ /2 ©3 r(* 2 ,i ■ 

2,i 2,i 


1 / 2 ) 


1 + 


*2.i®2,i 


-(*2,1 + 1) 


(13) 


We assume that quasineutrality prevails outside of perturbed De¬ 
bye regions, i.e., n 2fi = n 2p . 

Any metallic body embodied in the heliosheath plasma 
downstream of the shock charges up to an electric equilibrium 
potential <t >2 , which guarantees equal fluxes of ions and elec¬ 
trons reaching the metallic surface of this body per unit time 
(geometry taken to be planar). This behavior leads to the fol¬ 
lowing requirement (after dropping downstream indices “2” for 
simplification): 


■ Iff' 


j3 e (<t>) I I I (vcos6O/ e (v)v 2 dvd(/sin0d0 




(v cos 6)f{v)v~ dv dip sin 6 dO, (14) 


where /T( ( l ) ) and /?i(4>) denote the Boltzmann screening factors 
for electrons and ions, respectively. These factors describe the 
fraction of particles that can reach the wall against the electric 
potential <I>. The assumption of isotropic distribution functions 
leads to 


•f 


Pem /e(v)v 3 dv=j8 i«D) /i(v)v 3 dv. 




(15) 


We expect that the resulting equilibrium potential C> only affects 
the lowest-energy part of the distribution functions. Therefore, 
the Gaussian core of the kappa-distributed particles is the only 
screened population, leading to 


y8 e (0) _ exp (+2ech/m® 2 ) 
M®) exp (-2 e<D/M©f) 
which can be rewritten as 
1 1 3 


exp 


2 ed> 


, m@ 2 M® 2 


_ f f(v)v 3 dv 
f / e (v)v 3 dl-’ 

ff(v)v 3 dv 

f /e(v)v 3 dv 


(16) 


(17) 
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T e . Maxwellian thermal electrons lead to the classical Debye¬ 
screening length of /Id = yjKT e /4nne 2 . The effect of degen¬ 
erating Debye len gths has been recognized and emphasized by 
iTreumann et al.l (120041) . finding that the resulting Debye length 
dp may easily increase by factors of /I^/Td > 10 6 for low elec¬ 
tron kappa-indices of K e - 1.5. With a more relaxed approach, 
yet along the lines of these authors’ discussion, we obtain the 
following very similar conclusions. 

For the description of the effective screening of kappa elec¬ 
trons, we simply replace the Maxwellian temperature T e by the 
corresponding electron kappa temperature T£, which is given by 


KTl 


»e 


m Q 2 K e 
2 e K e - 3/2' 


( 22 ) 


Consequently, we find 


Fig. 2. Equilibrium potential of a metallic body in the presence of 
kappa-distributed ions and electrons according to Eq. CD- We show 
the potential as a function of K 2 ,c for different values of Af 2 ,i - We assume 
Ui = 400 km/s. 


We solve the remaining integrals in the above expression with 
02 .i/02.e ~ m/M, leading to 

f /j(v)v 3 dv rfa + 1) FO- e - 1/2) ytJTi (ire - 1) 

//e(t ; )v 3 dv r(/Ci - 1/2) r(K e + 1) y/ Ki M(Ki - 1) 

We find (with fl = 0; e /© 3 e ~ 1) 


__ / TV] j( 

D = y[r*JT e = \-®l—±—IKT e . (23) 

Assuming that the core of the kappa distribution is identical to 
the Maxwellian core (i.e., m&l = 2KT e ), we obtain the follow¬ 
ing result for the effective Debye length: 


dn = T D 


3/2 


(24) 


This modified Debye length has an interesting effect on the 
propagation of plasma waves. The general dispersion relation for 
electron-acoustic plasma waves (e.g., see lChenl 19741 Eq. (4-48)) 
is given by 


exp 


2e<t> U\ 
mU\ I© 3 


mU \' 
M® 2 , 



2e ° 7 2 , *2\ 

exp 



(19) 


where and /./*. denote the upstream solar-wind electron and 
pick-up ion Mach numbers. These numbers are given by values 
of the order /ii e - 10 and p * ; - 3. Therefore, 


exp 


' 2e$> 
MU\ 



rui +1) r(^e -1/2) ^/KJ^l(K e -1) 

r(/fi - 1/2) F(/c e + 1) y[i<\M 0q - 1) 


( 20 ) 


which leads to the following potential 

MU\ T(/f, + 1) r(«b - 1/2) y/kJTt (k c - 1) 
“ 200e T(if, - 1/2) rU e + 1) ( K[ - 1) 


( 21 ) 


As a consistency check, we note that this expression leads to 
the classical plasma-physics formula for = (D c in the limit 
of Maxwellian distributions (i.e., K e = k x —> oo) with identical 
temperatures T[ = T e = MU^/200K. 

In Figure [2] we show the resulting electric potential <f> as 
a function of the prevailing kappa-index ^.e of the shock-heated 
downstream electrons. This profile shows that the expected equi¬ 
librium potential drops to values of O < -30 V in the range of 
expected indices 1.5 < K e <2. This potential does not allow 
electrons with energies below 30 eV to reach the detector. 


4. Degenerated Debye length 

The degeneration of the Debye length is a direct consequence of 
highly nonthermal kappa-type distribution functions. The elec¬ 
tric screening by plasma-electron distributions with a largely ex¬ 
tended power-law tail is significantly less efficient compared to 
the screening by Maxwellian thermal electrons with temperature 


U) 

k 




AT, 


(25) 


In the general case of 7/ » 7) , and for very small wavevector 
values of k = 2n/A <s 2n/A k n , this dispersion relation allows 
for a branch of electron-acoustic waves that propagate with a 
phase/group velocity of 


u> du> I KT£ I m 2 

~k ~ ~dk~ ~ ^2M &e K e -3/2' 


(26) 


This branch is a special kappa-mode propagating with a typi¬ 
cal phase or group velocity that directly depends on the elec¬ 
tron kappa index K e . Testing plasma acoustic waves in this range 
of large wavelengths should hence directly reveal the prevailing 
kappa index K e and, therefore, the character of the suprathermal 
downstream electrons. 

On the other hand, the limit of larger wavevectors with k = 
2n/A > 2 tx! A h V) allows for a branch of nonpropagating waves (i.e., 
standing oscillations) with plasma eigenfrequencies given by 


CO 

ATg 1 


A 

M H 

r 

With Eqs. (l22l) and (l24l>. we can write 


KT K e 

2 ^- 3/2 / Anne 2 


MAg ~ \ 

v M 


(27) 


(28) 


for these oscillations. Under these conditions, the plasma, sur¬ 
prisingly enough, does not oscillate with the electron plasma 
frequency u > e , but it does oscillate with the ion plasma frequency 


Article number, page 4 of\5\ 





























































Fahr et al.: Electrons downstream of the termination shock 


u) p , which is a phenomenon that only quite rarely occurs in na¬ 
ture. The plasma oscillations recently registered by Voyager-1 
iGurnett et alj|2013l) . which were interpreted as electron plasma 
oscillations, may perhaps be reinterpreted as this type of ion 
plasma oscillations. As such, they would allow us to infer en¬ 
vironmental plasma densities of the order n - (m/M) 0.1 crrT 3 < 
10~ 4 cm -3 . 


5. Calculation of n in view of the downstream 
electron instabilities 

In Section 0 we introduced the quantity IT (see Eq. ©), which 
denotes the ratio of the thermal core widths, II = 0; /©, . Until 
this point, we did not determine a reasonable value for II. Pre¬ 
vious expressions based on an upstream-downstream transfor¬ 
mation of thermal-core velocities dFahr & Siewertfe013 f) appear 
to be irrelevant since their calculation relies on the assumption 
that upstream core electrons are independently transformed sim¬ 
ply into downstream core electrons according to the Liouville- 
Vlasov theorem. In reality, however, all upstream electrons over¬ 
shoot to the downstream side where the action of instabilities, 
such as the two-stream instability or the Buneman instability, 
redistribute and isotropize them. In the case of the two-stream 
instability (e.g., see lChenll974l) . electrons can excite ion oscilla¬ 
tions as long as their velocities are greater than the thermal core 
velocities of the protons. This fast relaxation of the electron dis¬ 
tribution function toward the downstream ion distribution func¬ 
tion then leads to a quasiequilibrium distribution. 

Consistent forms of such quasi-equilibria between parti¬ 
cle distribution functions and turbulence power spectra have 
been investigated in Section 2.5 in iFahr & Fichtnerl (20111) 
and, for steady state conditions, by lYoonl (20 111 i2012h and 
IZaheer & Yoonld2013l) . In all of these cases, the asymptotic state 
results in kappa distributions. Also, in our case, as a result of 
a shock-induced electron injection with velocity-space diffusion 
and relaxation described by a phase-space transport equation of 
the type 


?£-L±hn d/e\ . /e-/r 
dt v 2 <5v \ " dv) r eD 


(29) 


we expect solutions in form of k appa distributions. In fact, as 
shown by Treumann et al.l d2004l) . this kind of transport equa¬ 
tions leads to a quasi-equilibrium distribution in the form of a 
kappa distribution with a thermal core given by the downstream 
ion velocities: 0 2 „ (f p ) — 2 / , 2 ,p/« 2 ,e w i- We finally find (using 
Eq. (4) in iFahr & Siewertlfeo 1 3l) 

r. r. 2PT„U\m 2 

n = 0? e /0 1 e = —- - = - [2A(S, a) + B(s, a)] (30) 

jP i, v n 2 ,ejn 9 


withA(i, a) = Vcos 2 a + s 2 sin 2 a and Bis, a) = s 2 /A 2 (s,a). 

The above expressio n for a perpendicular shock with s = 2.5 
(Rich ardson et alJl2008h leads to Ifia = zr/2) = 1.33. This fi¬ 
nally shows that the quantity II is, in fact, of order unity, veri¬ 
fying a posteriori all of our above results that we calculated for 
II = 1. Using this more precise value for II, Eq. ( 1 1 011 leads to a 
marginally different value for the kappa-index of & 2 ,e = 1.522. 


6. Summary and conclusions 

Previous studies suggest that strongly-heated solar-wind elec¬ 
trons should appear in measurements as accelerated suprather- 
mal particles downstream of the termination shock. However, 


these electrons were not observed by Voyager in the heliosheath. 
We investigate this apparent contradiction and find that he¬ 
liosheath electrons are distributed according to a kappa-type dis¬ 
tribution function with an extended suprathermal tail (see Figure 
1). These highly suprathermal kappa-distributed electrons lead 
to a strong negative charging of all metallic bodies exposed to 
this plasma environment, consequently also charging up the Voy¬ 
ager spacecraft. 

A spacecraft potential of the order -30 V, as calculated in 
Figure [2j has a significant effect on the Voyager electron mea¬ 
surements in the heliosheath. Under these conditions, it repels 
thermal electrons in the energy range below 30 eV, leading to 
an increase in the prev iously determined upper limit of 3 eV 
(iRichardson et al.l l2008) for the electron temperature. The ions, 
on the other hand, are accelerated into the Faraday cups. The dif¬ 
ference in the derived bulk speeds, however, is negligible: cor¬ 
rected for a -30 V potential, the radial downstream proton bulk 
velocity increases from 130 km/s to about 132 km/s. 

In addition, this suprathermal distribution of downstream 
electrons also results in an unusually enlarged Debye length. As 
a consequence of this effect, the phase velocity = oj/k of elec¬ 
trostatic plasma waves depends on the effective kappa tempera¬ 
ture of the electrons in the heliosheath plasma environment. The 
detection of these plasma waves allows us to infer the effective 
kappa electron temperature as an observable quantity. These dis¬ 
tributions also permit a type of nonpropagating standing waves 
with the ion plasma frequency cj p as their eigen frequency. 
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